(N 

o 



Shape derivatives of boundary integral oper- 
ators in electromagnetic scattering. Part II: 
^ ; Application to scattering by a homogeneous 
dielectric obstacle 



C3 



Martin Costabel and Frederique Le Louer 



Abstract. We develop the shape derivative analysis of solutions to the 
, problem of scattering of time-harmonic electromagnetic waves by a pen- 

— ■ etrable bounded obstacle. Since boundary integral equations are a clas- 

sical tool to solve electromagnetic scattering problems, we study the 
CNj ■ shape differentiability properties of the standard electromagnetic bound- 

ary integral operators. The latter are typically bounded on the space 
f^**) , of tangential vector fields of mixed regularity TH~ 2 (divr, T). Using 

Helmholtz decomposition, we can base their analysis on the study of 
pseudo-differential integral operators in standard Sobolev spaces, but 
we then have to study the Gateaux differentiability of surface differ- 
ential operators. We prove that the electromagnetic boundary integral 
operators are infinitely differentiable without loss of regularity. We also 
give a characterization of the first shape derivative of the solution of 
the dielectric scattering problem as a solution of a new electromagnetic 
scattering problem. 

Keywords. Maxwell's equations, boundary integral operators, surface 
differential operators, shape derivatives, Helmholtz decomposition. 



1. Introduction 

Consider the scattering of time-harmonic electromagnetic waves by a bounded 
obstacle Q in M 3 with a smooth and simply connected boundary T filled with 
an homogeneous dielectric material. This problem is described by the sys- 
tem of Maxwell's equations with piccewisc constant electric permittivity and 
magnetic permeability, valid in the sense of distributions, which implies two 
transmission conditions on the boundary of the obstacle guaranteeing the 
continuity of the tangential components of the electric and magnetic fields 
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across the interface. The transmission problem is completed by the Silver- 
Miiller radiation condition at infinity (see [27J and [H]). Boundary integral 
equations are an efficient method to solve such problems for low and high 
frequencies. The dielectric scattering problem is usually reduced to a system 
of two boundary integral equations for two unknown tangential vector fields 
on the interface (see [5] and [H])- We refer to [5] for methods developed by 
the authors to solve this problem using a single boundary integral equation. 

Optimal shape design with a goal function involving the modulus of 
the far field pattern of the dielectric scattering problem has important appli- 
cations, such as antenna design for telecommunication systems and radars. 
The analysis of shape optimization methods is based on the analysis of the 
dependency of the solution on the shape of the dielectric scatterer, and a 
local analysis involves the study of derivatives with respect to the shape. An 
explicit form of the shape derivatives is desirable in view of their implementa- 
tion in shape optimization algorithms such as gradient methods or Newton's 
method. 

In this paper, we present a complete analysis of the shape differentia- 
bility of the solution of the dielectric scattering problem and of its far field 
pattern, using integral representations. Even if numerous works exist on the 
calculus of shape derivatives of various shape functional [TTJ E2 E3 [2D > 
in the framework of boundary integral equations the scientific literature is 
not extensive. However, one can cite the papers [55], [21] and [20], where 
R. Potthast has considered the question, starting with his PhD thesis [32] . 
for the Hclmholtz equation with Dirichlct or Neumann boundary conditions 
and the perfect conductor problem, in spaces of continuous and Holder con- 
tinuous functions. Using the integral representation of the solution, one is 
lead to study the Gateaux differentiability of boundary integral operators 
and potential operators with weakly singular and hypcrsingular kernels. 

The natural space of distributions (energy space) which occurs in the 
electromagnetic potential theory is TH~ 3 (divr, T), the set of tangential vec- 
tor fields whose components are in the Sobolev space H~z(T) and whose 
surface divergence is in H-2(T). We face two main difficulties: On one hand, 
the solution of the scattering problem is given in terms of products of bound- 
ary integral operators and their inverses. In order to be able to construct 
shape derivatives of such products, it is not sufficient to find shape deriva- 
tives of the boundary integral operators, but it is imperative to prove that 
the derivatives are bounded operators between the same spaces as the bound- 
ary integral operators themselves. On the other hand, the very definition 
of shape differentiability of operators defined on the shape-dependent space 
TH 3 (divr, T) poses non-trivial problems. Our strategy consists in using the 
Hclmholtz decomposition of this Hilbert space which gives a representation 
of a tangential vector field in terms of (tangential derivatives of) two scalar 
potentials. In this way, we split the analysis into two steps: First the Gateaux 
differentiability analysis of scalar boundary integral operators and potential 
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operators with strongly and weakly singular kernels, and second the study of 
shape derivatives of surface differential operators. 

This work contains results from the thesis [21] where this analysis has 
been used to develop a shape optimization algorithm of dielectric lenses in 
order to obtain a prescribed radiation pattern. 

This is the second of two papers on shape derivatives of boundary in- 
tegral operators, the first one [H] being aimed at a general theory of shape 
derivatives of singular integral operators appearing in boundary integral equa- 
tion methods. 

The paper is organized as follows: 

In Section [2] we recall some standard results about trace mappings and 
regularity properties of the boundary integral operators in electromagnctism. 
In Section [3] we define the scattering problem for time-harmonic electromag- 
netic waves at a dielectric interface. We then give an integral representation 
of the solution — and of the quantity of interest, namely the far field of the 
dielectric scattering problem — following the single source integral equation 
method developed in jS]. 

The remaining parts of the paper are dedicated to the shape differen- 
tiability analysis of the solution of the dielectric scattering problem. We use 
the results of our first paper [S] on the Gateaux differentiability of boundary 
integral operators with pseudo-homogeneous kernels. We refer to this paper 
for a discussion of the notion of Gateaux derivatives in Frechct spaces and of 
some of their basic properties. In Section |4] we discuss the difficulties posed 
by the shape dependency of the function space TH _5 (divr,r) on which the 
integral operators arc defined, and we present a strategy for dealing with this 
difficulty, namely using the well-known tool [TO] of Hclmholtz decomposition. 
In our approach, we map the variable spaces TH _ 3 (divr r , IV) to a fixed ref- 
erence space with a transformation that preserves the Hodge structure. This 
technique involves the analysis of surface differential operators that have to 
be considered in suitable Sobolev spaces. Therefore in Section ([5]) we recall 
and extend the results on the differentiability properties of surface differen- 
tial operators established in Section 5]. Using the rules on derivatives of 
composite and inverse functions, we obtain in Section® the shape differentia- 
bility properties of the solution of the scattering problem. More precisely, we 
prove that the boundary integral operators are infinitely Gateaux differen- 
tiablc without loss of regularity, whereas previous results allowed such a loss 
|30j . and we prove that the shape derivatives of the potentials are smooth 
away from the boundary but they lose regularity in the neighborhood of the 
boundary. This implies that the far field is infinitely Gateaux diffcrcntiable, 
whereas the shape derivatives of the solution of the scattering problem lose 
regularity. 

These new results generalize existing results: In the acoustic case, using 
a variational formulation, a characterization of the first Gateaux derivative 
was given by A. Kirsch [5T] for the Dirichlet problem and then by Hettlich 
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|16H17j for the impedance problem and the transmission problem. An alterna- 
tive technique was introduced by Kress and Paivarinta in [53] to investigate 
Frechet differentiability in acoustic scattering by the use of a factorization of 
the difference of the far-field pattern of the scattered wave for two different 
obstacles. In the electromagnetic case, Potthast used the integral equation 
method to obtain a characterization of the first shape derivative of the solu- 
tion of the perfect conductor scattering problem. In [22) , Kress improved this 
result by using a far-held identity and in [H] Kress and Haddar extended 
this technique to acoustic and electromagnetic impedance boundary value 
problems. 

At the end of Section H] we obtain a characterization of the first shape 
derivative of the solution of the dielectric scattering problem as the solu- 
tion of a new electromagnetic transmission problem. We show by deriving 
the integral representation of the solution that the first derivative satisfies 
the homogeneous Maxwell equations, and by directly deriving the boundary 
values of the solution itself we see that the first derivative satisfies two new 
transmission conditions on the boundary. 

In the end we will have obtained two different algorithms for computing 
the shape derivative of the solution of the dielectric scattering problem and of 
the far field pattern: A first one by differentiating the integral representations 
and a second one by solving the new transmission problem associated with 
the first derivative. 

The characterization of the derivatives as solutions to boundary value 
problems has been obtained in the acoustic case by Kress [6] , Kirsch [21] , Hct- 
tlich and Rundell [T8] and Hohage [T^] and has been used for the construction 
of Newton-type or second degree iterative methods in acoustic inverse obsta- 
cle scattering. Whereas the use of these characterizations requires high order 
regularity assumption for the boundary, we expect that the differentiation of 
the boundary integral operators does not require much regularity. Although 
in this paper we treat the case of a smooth boundary, in the last section we 
give some ideas on possible extensions of the results of this paper to non- 
smooth domains. 

A final remark on the terminology of derivatives with respect to the 
variation of the domain: For the solutions of boundary value problems with 
a moving boundary, in the context of continuum mechanics one distinguishes 
frequently between material derivatives and shape derivatives. The former 
correspond to the Gateaux derivative with respect to the deformation, which 
is then interpreted as a flow associated with a velocity field, when the so- 
lution is pulled back to a fixed undeformed reference domain (Lagrangian 
coordinates). Depending on the interpretation of the flow, other names are 
used, such as "Lagrangian derivative" , or "substantial derivative" . The shape 
derivatives, on the other hand, correspond to the Gateaux derivative of the 
solution in the deformed domain (Eulerian coordinates) . One also talks about 
"Eulerian derivative" . The difference between the two has the form of a con- 
vection term, which can lead to a loss of one order of regularity of the shape 
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derivative on the support of the velocity field [25] . In the context of this ter- 
minology, the derivatives of the boundary integral operators that we studied 
in Part I of this paper [3] for general pseudohomogeneous kernels and in this 
paper for the boundary integral operators of electromagnetism, correspond 
to the Lagrangian point of view, because they are obtained by pull-back to 
a fixed reference boundary. However, the derivative of the solution of di- 
electric scattering problem that we construct in Section - with the help 
of the derivatives of the boundary integral operators - is the Eulerian shape 
derivative. We also observe the loss of one order of regularity of the shape de- 
rivative with respect to the solution of the transmission problem. Finally, for 
the far field, the two notions of derivative coincide, because the deformation 
has compact support. 

2. Boundary integral operators and their main properties 

Let O be a bounded domain in R 3 and let ft c denote the exterior domain 
]R 3 \f2. Throughout this paper, we will for simplicity assume that the bound- 
ary r of is a smooth and simply connected closed surface, so that ft is 
diffcomorphic to a ball. 

We use standard notation for surface differential operators and bound- 
ary traces. More details can be found in For a vector function v E 
^(M 3 , C 3 ) with k e N*, we denote by [Vu] the matrix the i-th column of 
which is the gradient of the i-th component of v, and we set [Dv] = [Vt>] T . 
Let n denote the outer unit normal vector on the boundary T. The tangential 
gradient of a complex- valued scalar function u € < ^ fc (r,C) is defined by 

V r u = V{t|r - (Vii|r • n) n, (2.1) 

and the tangential vector curl is defined by 

curlr u = Vtt|r A n, (2-2) 

where u is a smooth extension of u to the whole space R 3 . For a complex- 
valued vector function u G ^ fc (r,C 3 ), we denote [Vrw] the matrix the i-th 
column of which is the tangential gradient of the i-th component of u and 
we set [Dr«] = [Vrit] T - 

The surface divergence of u £ ^ fc (r,C 3 ) is defined by 

divr it = divii|r — ([Vit|r]'n. • w) , (2-3) 

and the surface scalar curl is defined by 

curlr u — n ■ (curl u) . 

These definitions do not depend on the choice of the extension ii. 

Definition 2.1. For a vector function v G c rf ca (^l,C 3 ) a scalar function v E 
^°°(f2, C) and k E C \ {0}, we define the traces : 
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InV = T^V = n ■ (VU)| r , 

-fr>v := n A (v)\ T (Dirichlet), 
7at„ , u := — n A (curl 17)1 (Neumann). 

K 

We define in the same way the exterior traces 7 C , 7°, -f D and 7^ . 

For a domain G C K 3 we denote by H S {G) the usual L 2 -based Sobolev 
space of order s g R, and by H? oc (G) the space of functions whose restric- 
tions to any bounded subdomain B of G belong to H S (B). Spaces of vector 
functions will be denoted by boldface letters, thus 

H S (G) = {H s (G)f . 

If D is a differential operator, we write: 

H s (D,fi) = {«eH s (fl):DH€H s (!l)} 

Hf oc (D,IF) = {it e Hf oc (IF) : Du G Hf oc (TF)}. 

The space H s (D,f2) is endowed with the natural graph norm. When s = 0, 
this defines in particular the Hilbert spaces H (curl, fl) and H (curl curl, ft). 

We introduce the Hilbert spaces H S (T) = j^H s+ ^(Q)^ , and TH s (r) = 
7d ^H i+ 2(f2)^. For s > 0, the trace mappings 

7 : H s+ i{n) ^ H S (T), 

lD : H s+ 5(fi) TH s (r) 

are then continuous. The dual of H S (T) and TH s (r) with respect to the L 2 
(or L 2 ) scalar product is denoted by H~ S (T) and TH~ s (r), respectively. 

The surface differential operators defined above can be extended to 
Sobolev spaces: For s€l the tangential gradient and the tangential vector 
curl are obviously linear and continuous operators from H S+1 (T) to TH s (r). 
The surface divergence and the surface scalar curl can then be defined on tan- 
gential vector fields by duality, extending duality relations valid for smooth 
functions 

J (div r j) • fds = - J j ■ \7 r ipds for all j e TH s+1 (r), <p e £T s (r), 

(2.4) 

J (curlrj) • <pds = J j ■ curl r <pds for all j e TH s+1 (r), ip € ff~ s (r). 

(2.5) 

We have the following equalities: 

curlr Vr = and divr curlr = (2-6) 
divr(n A j) = — curlr j and curlr (n A j) = divr j (2.7) 
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Definition 2.2. Let s £ R. We define the Hilbert space 

TH s (div r ,r) = {j e TH s (r) ; div r j e tf s (r)} 

endowed with the norm 

II • llTH s (div r ,r) = (ll ' llTH*(r) + 1 1 div r • I \h> (r)) ■ 

Lemma 2.3. The operators 7d and are linear continuous from ( ^' 00 (f2, C 3 ) 
to TL 2 (r) and they can be extended to continuous linear operators from 
H(curl,f2) and H(curl, f2)n H(curl curl, f2), respectively, to TH 5 (divr, L). 

For u £ Hi oc (curl, ft c ) and v £ Hi oc (curl curl, f2 c )) we define j^u and 
in the same way and the same mapping properties hold true. 
Recall that we assume that the boundary F is smooth and topologically 
trivial. For a proof of the following result, we refer to [T| I7j [28] . 

Lemma 2.4. Let s £ R. The Laplace- Beltrami operator defined by 

Apw = divr Vr« = — curlr curlr u. (2-8) 

is linear and continuous from H S+2 (T) to H S (T). For f £ H S (T) and u £ 
H S+2 (T), the equation Apu = f has the equivalent formulation 

J VrM-Vr^ds = - J f ■ tpds, for all ip £ H- S (T). (2.9) 

The operator Ay ■ H S+2 (T) — > H S (T) is Fredholm of index zero, its kernel and 
cokernel consisting of constant functions, so that Ar '■ H S+2 (T)/C — > H%(T) 
is an isomorphism. Here we define the space H%(T) by 

f £ H*(T) <=► f£ H S (T) and J^fds = 0. 

For f £ H'(T) we denote the unique solution u £ H S+2 (T)/C of (|2.9|l by 
u = Af 1 f. 

This result is due to the injectivity of the operator Vr from H S+2 (T)/ 
C to TH s+1 (r), the Lax-Mil gram lemma applied to (I2.9[) for s — — 1, and 
standard elliptic regularity theory. Note that curlr is also injective from 
H S+2 (T)/C to TH s+1 (r), and by duality both divr and curlr are surjective 
from TH s+1 (r) to H*(T) 

Notice that curlr is defined in a natural way on all of H s+1 (r) and 
maps to H%(T), because we have curlr (ipn) = for any scalar function <p £ 
H S+1 (T). Thus ()2.5p is still valid for a not necessarily tangential density 
j £ H s+ (r). An analogous property for divr defined by (|2.3p is not available. 

We now recall some well known results about electromagnetic potentials. 
Details can be found in [TJ E9 U |2D1 121] • 

e in\x-y\ 

Let k be a positive real number and let Gain, \x — y\) = — ; r 

4ir\x — y\ 

be the fundamental solution of the Helmholtz equation Au + n 2 u = 0. The 
single layer potential ^ K is given by 
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4Vu(z) = J^G a {K,\x-y\)u{y)ds{y) x eR 3 \T, 
and its trace by 

V K u(x) = l G a (n,\x - y\)u(y)ds(y) xeT. 



As discussed in the first part of this paper [5], the fundamental solution is 
pseudo- homogeneous of class — 1. The single layer potential ^ K u is continuous 
across the boundary T. As a consequence we have the following result : 

Lemma 2.5. Let s€l. The operators 

* K : H s -i(T) -> ff^fi) n Kf(SF) (ff-*(T) -> i/ 1 s + 1 (K 3 ) « < ±) 

F« : H s -i(T) -> H s +5(r) 
are continuous. 

The electric potential operator 4^ is defined for j 6 TH 5 (divr, T) 

by 

*i3„j := K& K j + -V* K div r j. 
In R 3 \ r, this can be written as $e j := — curl curl because of the 

K 

Helmholtz equation and the identity curl curl = — A + Vdiv. 

The magnetic potential operator ^m k is defined for m € TH 5 (divr, T) 

by 

^>M K 'm '■= curl^Km. 

We denote the identity operator by I. 

Lemma 2.6. The potentials operators and ^m k are continuous from 

TH - '(div r ,r) to Hi oc (curl, M 3 ). For j £ THT^ (div r , T) we have 

(curl curl -K 2 l)^ E J = and (curl curl -k 2 !)^ Mri rn = in R 3 \r, 
and ^E K j and ^ m^tti satisfy the Silver-Miiller condition. 

We define the electric and the magnetic far field operators for a density 
j and an element x of the unit sphere S 2 of R 3 by 



(2.10) 



*H Ox) = K & A / e-^ x -yj{y)ds{y) A x 



r 



j(x) = iKxA[ / e-^j(y)ds(y) 



r 



These operators are bounded from TH s (divr,r) to 
TL 2 ^ 2 ) ={/i£ L 2 (S* 2 ); • 2 = 0}, for all s e R. 

We can now define the main boundary integral operators: 

C K j{x) = - J n(x) A curlcurl a; {G' Q (K, \x ~ y\)j(y)}ds(y) 

-K nAV K j + -curl r V fi div r j)(x), (2.11) 

K J 
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and 



with 



n(x) A curl a {G (K, \x - y\)j(y)}ds(y) 



= (D K j-B K j)(x), 



B K j(x) = J V x G a (K, \x - y\) (j(y) ■ n{x)) ds(y), 
D K j{x)= [ (V*G a (K,\x-y\)-n(x))j(y)ds(y). 



(2.12) 



The operators M K and C K are bounded operators from TH 2 (divr,T) 
to itself. 



3. The dielectric scattering problem 

We consider the scattering of time-harmonic waves at a fixed frequency lu 
by a three-dimensional bounded and non-conducting homogeneous dielectric 
obstacle represented by the domain Q. The electric permittivity e and the 
magnetic permeability \x are assumed to take constant positive real values 
in fi and tt c . Thus they will be discontinuous across the interface T, in gen- 
eral. The wave number is given by k — lj^JJU. We distinguish the dielectric 
quantities related to the interior domain £1 through the index i and to the 
exterior domain Q c through the index e. The time-harmonic Maxwell sys- 
tem can be reduced to second order equations for the electric field only. The 
time-harmonic dielectric scattering problem is then formulated as follows. 

The solution of the dielectric scattering problem : Consider the scat- 
tering of a given incident electric wave E mc 6 Hi oc (curl, R 3 ) that satisfies 
curl curl E mc — K^E mc = in a neighborhood of tt. The interior electric field 
E* £ H(curl, tt) and the exterior electric scattered field E s G Hi oc (curl, Q c ) 
satisfy the time-harmonic Maxwell equations 



2£S 



K 



curl curl E 
curl curl E l 

the two transmission conditions, 

n A E' = n A (E 
[i l ~{n A curl E l ) = fj,~ 
and the Silvcr-Miiller radiation condition: 



= 
= 

1 (n A curl E*) = A curl(E s 



lim | a; | 



curl E s (x) A 



in O, 

in n c 



on r 
on r 



0. 



(3.1) 
(3.2) 

(3.3) 
(3.4) 

(3.5) 



It is well known that the problem (|3.1j) - (|3.5j) admits a unique solution 
for any positive real values of the exterior wave number n e . We refer the 
reader to [5J |H1 for a proof via boundary integral equation methods. 
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To analyze the dependency of the solution on the shape of the scatterer 
SI, we will use an integral representation of the solution, obtained by the 
single boundary integral equation method developped by the autors in [5] . It 
is based on the layer ansatz for the exterior electric field E s : 

E S = -^e k J - iV*M Ke C^j in M 3 \ H (3.6) 

where 77 is a positive real number, j 6 TH~ 5 (divr,r) and the operator Cq 
is defined for j € TH 5 (div r , T) by 

CqJ = -nAV j - curl r V div r j. 

Thanks to the transmission conditions and the Stratton-Chu formula, we 
have the integral representation of the interior field 

E* = --{^ E Al C N^ mC + N e j}) - (* Mk . {7dE" c + L e j}) in SI (3.7) 
p 

Kip e 

where p = and 

K e pi 

L e = C lt .+iri(-~I + M K }(%, 

N e = (-^1 + M K )j + ir]C Ke CZ. 

The exterior Dirichlet trace applied to the right-hand side Q3.7j) vanishes. The 
density j then solves the following boundary integral equation 



Sj = p (-il + M Ki ) L e j+C Ki N e j = -p (-±1 + M Kz ) lD E mc + C RtlN ^ E mc . 

(3.8) 

Theorem 3.1. The operator S from TH~ 5 (divr, T) to itself is linear, bounded 
and invertible. Moreover, given the electric incident field E mc e Hi oc (curl, R 3 ) ; 
the integral representations (|3.6|) . (|3.7I) of E l and E s give the unique solution 
of the dielectric scattering problem for all positive real values of the dielectric 
constants pi, p e , e; and e e . 

An important quantity that is of interest in many shape optimization 
problems is the far field pattern of the electric field, defined on the unit sphere 
of R 3 by 




We have E°° G TL 2 (5 2 ) n ^ °°(S 2 , C 3 ). To obtain the integral representation 
of the far field E°° of the solution, it suffices to replace in ()3.6|) the potential 
operators and ^m k by the far field operators Vt'^ 3 and ^fj defined 

in (|2.10j) . respectively. 

In the method we have described, the solution E = (E 4 , E s ) and the far 
field E°° are constructed from operators defined by integrals on the boundary 
F and the incident field. For a fixed incident field and fixed constants Hi, n e , 
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/Xj, /i e , these quantities therefore depend on the geometry of the boundary 
T of the scatterer f2 only. In the sequel we analyze the T-dependence of the 
solution following the definition of shape derivatives and the notations of 
section 4 of the paper [9] . 

4. Shape dependence via Helmholtz decomposition 

Let us fix a reference domain Q. We consider variations of Q generated by 
transformations of the form x i— > x+r(x) of points x € R 3 , where r is a smooth 
vector function defined on T. This transformation deforms the domain fi in 
a domain f2 r of boundary T r . The functions r are assumed to belong to the 
Frechet space <^°°(r,R 3 ). For e > and some metric on ^f°°(r,R 3 ), we 
set 

B°°(0,e) = {re ^°°(r,R 3 ), ^(O.r) < e} . 

In the following, we choose e small enough so that for any r <E B°°(0,e), 
(I + r) is a diffcomorphism from r to T r = (l + r)T = {x r = x + r(x)\ x € F}. 

The aim of this paper is to study the shape differentiability, that is, the 
Gateaux differentiability with respect to r, of the functionals mapping r to 
the solution 

(^(r),^(r)) = (E*(r P ) J E'(T P )) 

of the dielectric scattering problem with obstacle Vt r , and to the far field 
<?°°(r) = E°°(r r ). 

In the following, we use the superscript r for integral operators and trace 
mappings pertaining to T r , while functions defined on T r will often have a 
subscript r. According to the boundary integral equation method described 
in the previous section, we have 

S s (r) = (-*^ - i V W MKe C*' r ) jr ^ ^ r = M?\W r , (4.1) 

where j r solves the integral equation 

ST • ( ^ j \ n \ r rinc fir r rinc 

Jr = ~P ( ~2 I + M «i J TD E - C Ki lN Ke E . 

and 
with 

<g tot (r) = e nc + £"(r). (4.3) 

The far field pattern of the dielectric scattering problem by the interface F r 
is 

^°°M = (-*SI-«;:^ r )^- 

As defined in (|3.8|) . the operator S r is composed of the operators C T K , MJ1 , 
C£. et M£., which are all bounded operators on the space TH" (divr r ,r r ). 
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Therefore we have to study the Gateaux differentiability of the following 
mappings on B°°(0,e) 

r^M£,C£eJ^(TH-*(div rr ,r r )) 

y r M K ,* r En G J2?(TH-i(divr r ,r r ),H(curl,fi r )UH Ioc (ciirl,n«)) 
r i— > ^^' .'1'^;' G ^(TH-^(divr r ,r r ),TL 2 (5 2 ) n ^°°(S 2 ,C 3 )). 

Finally, the differentiability properties of the mapping r \— > Cq t can be de- 
duced from those of the mapping r i— > C r K . 

In this approach, several difficulties have to be overcome. The first one 
is that if we want to find the derivatives of the solution of the scattering 
problem, which is given as a product of operators and of their inverses, all 
defined on the same space TH~ 5 (divr, T) (for the derivative at r = 0), it 
is necessary to prove that the derivatives themselves are defined as bounded 
operators on the same space, too. On the other hand, the very definition 
of the differentiability of operators defined on TH~ 1 (divr, T) raises non- 
trivial questions. For reducing the variable space TH~ 5 (divr r , T r ) to a fixed 
reference space, it is not sufficient, as we did in the scalar case studied in 
the first part [§], to use a change of variables. Let us discuss this question: 
How to define the shape derivative of operators defined on the variable space 
TH"^(div rr .,r r )? in detail. 

We recall the notation r r for the "pullback" induced by the change of 
variables. It maps a function u r defined on r r to the function r r u r = u r o(l+r) 
defined on V. For r G B°°(0, e), the transformation r r is an isomorphism from 
F*(r r ) to H\T). We have 

(r r u r )(x) = u r (x + r(x)) and (r~ 1 u)(a; r ) = u(x). 

The natural idea to use this for a product of operators, proposed by Potthast 
in [22] in the acoustic case, is to insert the identity T~ Y T r = i ^ between 

the factors. This allows to consider integral operators on the fixed boundary 
r only and to would require study the differentiability of the mappings 

r^r r C r K T-\ r^TrM^T- 1 , r^*^" 1 , (4.4) 

but as has been already pointed out in |30) . difficulties remain. The main 
cause for this is that r r does not map vector fields tangential to T r to vector 
fields tangential to T, and in particular, 

r r (TH-^ (divr,. , r r )) ? TH^ (div r , T). 

This will lead to a loss of regularity if we simply try to differentiate the 
mappings in (|4.4[) . Let us explain this for the operator M K . The operator M£ , 
when acting on vector fields tangential to r r , has additional regularity like 
what is known for the scalar double layer potential, namely it has a pseudo- 
homogeneous kernel of class —1, whereas it is of class when considered on 
all vector fields (see (|2.12p 'l. If we differentiate the kernel of Tr-M^r" 1 , we 
will not obtain a pseudo-homogeneous kernel of class —1 on the set of vector 
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fields tangential to F, so that we find a loss of regularity for the Gateaux 
derivative of T r M r K T~ x . 

For mapping tangent vector fields to tangent vector fields, the idea of 
Potthast was to use projectors from one tangent plane to the other. Let 
us denote by 7r(r) the pullback r r followed by orthogonal projection to the 
tangent plane to T. This maps any vector function on T r to a tangential 
vector function on T, and we have 

(7r(7")u r )(x) = u r (x + r{x)) — (n(x) ■ u r (x + r(x))) n(x). 

The restriction of 7r(r) to tangential functions on T r admits an inverse, de- 
noted by 7r (r), if r is sufficiently small. The mapping w (r) is defined 
by 

{-K-\r)u){x + r{x)) = u{x) - n(x) ^ + ^ ' , 
v w n K " w v ' n r (x + r(x)) ■ n(x) 

and it is easy to see that 7r(r) is an isomorphism between he space of con- 
tinuous tangential vector functions on r r and on T, and for any t between 
TH*(T r ) and TH*(T). 

In the framework of continuous tangential functions it suffices to insert 
the product 7r~ 1 (r)7r(r) = I-nf°(rv) between factors in the integral represen- 
tation of the solution to reduce the analysis to the study of boundary integral 
operators defined on T < ^ p °(r), which does not depend on r. In our case, we 
would obtain operators defined on the space 

7r(r) (ThT5(di Vrv ,rv)) = |« G TH~*(r),divr r (7r _1 (r)«) G ff"5(r r )}, 

which still depends on r and is, in general, different from TH~ 5 (divr, T). 

We propose a different approach, using the Helmholtz decomposition 
of the space TH~ 5 (divr r ,r r ) to introduce a new pullback operator P r that 
defines an isomorphism between TH" (divr r , T r ) and TH~ 5 (divr, T). 

Recall that we assume that the boundary T is smooth and simply con- 
nected. We have the following decomposition. We refer to [TU] for the proof. 

Theorem 4.1. The Hilbert space TH~ 5 (divr, T) admits the following Helmholtz 
decomposition: 

TH"5(div r ,r) = Vr-H"'(r) ©curl r if5(T). (4.5) 

Since e is chosen such that for all r £ B°°(Q,s) the surfaces T r are 
still regular and simply connected, the spaces TH~ 5 (divr r ,r r ) admit similar 
decompositions. 

The operator of change of variables r r is an isomorphism from Hi(T r ) 
to H 2 (r) and from H 2 (T r ) to H 2 (T), and it maps constant functions to con- 
stant functions. Let j r e TH _ 2 (divr r , T r ) and let j r = Vr r p r +curlr,. q r be 
its Helmholtz decomposition. The scalar functions p r and q r are determined 
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uniquely up to additive constants. The following operator : 

P r : TH"5(div rr ,r r ) — > TH _ *(div r ,r) 

j r = V r ,. p r + curl Fr q r i-> j = V r (r r p r ) + curl r {r r q r ) 

(4-6) 

is therefore well defined, linear, continuous and invertible. Its inverse P r 1 is 
given by 

P- 1 : TH-^(div r ,r) ->■ TH-5(divr, ., T r ) 

j = V r p + curlr q >-> j r = Vr r T~ x {p) + curlr, T~ x {q). 

(4.7) 

Obviously for r = we have P r = P_ 1 = I i . We can now insert 

J r TH~2(div r ,r) 

the identity I _i = P~ P r between factors in the integral repre- 

TH 2 (divr r ,r,.) 

sentation of the solution (^"(r), S s (r)), and we are finally led to study the 
Gateaux differentiability properties of the following mappings, defined on 
r-independent spaces. 

B°°(Q,e) jSf(TH-5(di Vr ,r),H(curl,ifp)) : r ^P' 1 

B°°(0,e) -> JSf(TH~*(divr,r),H(curl,Xp)) : r^^P^ 1 
B°°(Q,e p ) -> if(TH s (div r ,r),TH-5(div r ,r)) : r H- Pr-M^P" 1 
£°°(0,£ p ) -> ^(TH-5(divr,r),TH^5(div r ,r)) : r P^P" 1 

(4.8) 

where i4T p is a compact subset of M 3 \r. These mappings are composed of 
scalar singular integral operators, the shape derivatives of which we stud- 
ied in the first part [S], of surface differential operators, and of the inverse 
of the Laplace-Beltrami operator, which appears in the construction of the 
Hclmholtz decomposition. 

Let us look at the representation of the operators in (|4.8[) in terms of 
the Helmholtz decomposition. 
Helmholtz representation of ^ r E P~ 

The operator P r -1 is defined for j — Vr p + curl r q G TH~^(div r ,r) 

and x G by: 

W E P- X j(x) = kJ^ G a (K, \x - y r \) {V Tr .r- l p) (y r )ds(y r ) 

+ kJ G a (K, \x - y r \) (curlr,. T~ l q) {y r )ds{y r ) 
+ -V / G a (n, \x - y r \) (Ar r r~V) (Vr)ds(y r ). 

Helmholtz representation of ty r M P^ 1 

The operator P r _1 is defined for j = Vr p + curlr 9 G TH~^ (div r , T) 
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and x £ K p by: 

*M„ P r 1 3{x) = cur1 ^ G a (K, \x - y r \) {Vr^p) (y r )ds(y r ) 

+ curl J G a (n, \x - y r \) (curlrvT" 1 ^) (y r )(y r )ds(y r ). 

Helmholtz representation of P r C£P~ 

Recall that for j r e TH" 5 (divr r , r r ), the operator C£ is defined by 



Gr Jr{ x r)= -nn r (x r ) A I G a {n,\x r - y r \)j r (y r )ds(y r ) 
1 



rv 



n r (a; r )AVp r / G a [n, \x r - y r \) div rr j r (y r )ds(y r ). 
K r Jr r 

We want to write G T K j T in the form Vr r P r + curlr r Q r - Using formulas 
(pTH) . we find 

div Fr C r J r = A r ,,P r and curlr r C r J r = -A Fr Q r - 

As a consequence we have for x r G T r 

Pr(x r ) = -« Af^divr,. fn r (x r ) A J G a {K,\x r - y r \)j r (y r )ds(y r ) 

(4.9) 

and 



Qr(ir) = -K (-Ap^ 1 ) curlr, (n r (x r ) A J G a (K,\x r - y r \)j r (y r )ds(y r ) 
--(-Ar r ) curl r ,.(- curl r ,.) / G a (n,\x r - y r \) dw Fr j r (yr)ds(y r ) 
= k Ap 1 curlr r ^n r {x r ) A J G a (K,\x r - y r \)j r (y r )ds(y r )^j 
+ - I G a (K, \x r - y r \) div Fr j r (y r )ds(y r ). 

The operator P^P,: 1 is defined for j = V r p + curl r q € TH - *(div r ,r) 

by 

PrC^; 1 = V r P(r) + curlr Q(r), 

with 
P(r)(x) = 

-k (r r A^ divr r t^T 1 ) ^(r r n r )(a;)AT r | / <3 (k, | • -y r \)(Vr r T~ 1 p){y r )ds(y r ) 



+ Ga(n, | • -y r |)(curlr r r r 1 <j)(2/r)rfs( 2 / r )|(.T)^ 
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and 

Q(r)(x) = 

k (r r A^ curlr r r^ 1 ) ((T r n r )(x) A r r | J G a (n, \ ■ -y r \){^r r T^ 1 p){y r )ds{y r ) 
+ J G a (K, I • -y r \)(curlr r T^ 1 q)(y r )ds(y r )^(x)j 

+ ~ Tr {J G °-( K ' \ ' -yr\){&r r Tr 1 p){y r )ds{y r )J{x)- 

Helmholtz representation of PrM^P^ 1 

Recall that for all j r e TH^ (divr r ,r r ), the operator M r K is defined by 



M Ur( x r) = I ((V Xr G a (n, \x r - y r \)) ■ n r (x r ))j r (y r )ds(y r ) 
V Xr G a (K, \x r - y r \){n r (x r ) ■ j r (y r ))ds(y r ). 
Using the equalities (|2.7I) and the identity curl curl = — A + Vdiv, we have 
divrv M r J r {x r ) = n r {x r ) J curlcurr r {G a (n, \x r ~ y r \)j r (y r )} ds(y r ) 



= K 2 n r (x r ) ■ I {G a (K,\x r - y r \)j r (y r )}ds(y r ) 
8 



dn r {x r ) 



{G a (n, \x r - y r \) divrv j r (y r )} ds(y r ). 



Proceeding in the same way as with the operator P r C^P r , we obtain that 
the operator P r M£P~ is defined for j = Vr p + curlr q G TH 5 (divr, T) 
by: 

P r M r K P^j = V r P'(r) + curlr Q'(r), 



with 
P'(r)(x) = 



|k 2 J n r ■ {G a (K, I • -y r \) curlr,. r r 1 q(y r )} ds(y r 



k 2 J n r ■ {G a (K, | • -2/ r |)Vr r T r 1 p(y r )} ds(y r ) 

+ -^ G ^ I ' -yr\)(^r r T 7 7 1 p)(y r )ds(y r ) \ (.ri. 
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and 

(tvA^ 1 curlrv t^ 1 ) t t I^J (VG„(k, | • -y r \) ■ n r ) (curl r ,. T~ 1 q)(y r )ds{y r ) 
+ ((VG„(k, I • -y r \) ■ n r ) (Vr r r r - l p)(y r )ds(y r ) 
- J VG„(k, I • -y r \) (n r ■ (curlrv T-T^Xz/r)) ds{y r ) 

- jf VG„(k, I • -y r |) (n r ■ (Vr r r"V)(yr)) da(lfr) J (»)• 

These operators are composed of boundary integral operators with weakly 
singular kernel and of the surface differential operators defined in section 
[2J Each of these weakly singular boundary integral operators has a pseudo- 
homogeneous kernel of class -1. The ^"^-Gateaux differentiability properties 
of such boundary integral operators has been established in the preceding 
paper [S]. It remains now to show that the surface differential operators, 
more precisely rrVr^TjT 1 , r r curlr,. r" 1 , i~ r divr r r" 1 , r r curlr r t" 1 , as well 
as r r Ar r T~ 1 and its inverse, preserve their mapping properties by differenti- 
ation with respect to r. 



5. Gateaux differentiability of surface differential operators 

The analysis of the surface differential operators requires the differentiability 
properties of some auxiliary functions, such as the outer unit normal vector 
n r and the Jacobian J r of the change of variable x i— > x + r(x). We recall 
some results established in the first part j9] section 4]. For the definition of 
Gateaux derivatives and the corresponding analysis, see |33j . 

Lemma 5.1. The mapping Af : B°°(0,e) 3 r H> r r n r = n r o (I + r) € 
c ta°°(T,M 3 ) is 't? oa -Gateaux- differentiable and its first derivative in the direc- 
tion o/£ G <tf°°(r,M. 3 ) is given by 

dM[r,t} = -[r r Vr r (Tr 1 0]M(r). 

Lemma 5.2. The mapping J from r G S°°(Q,e) to the surface Jacobian 
J r G ^°°(r,R) is "tf 00 -Gateaux differentiable and its first derivative in the 
direction o/£ € < if 00 (r,IR 3 ) is given by 

dj[r a ,£} = J ro ■ (r ro divr^Cr" 1 ^))- 

The differentiability properties of the tangential gradient and of the 
surface divergence in the framework of classical Sobolev spaces is established 
in section 5]. 
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Lemma 5.3. The mapping 

Q: B°°(0, e) -> ^f(H s+1 (T), H s (r)) 
r i — ^ T r Vr r T r _1 

is ff 00 -Gateaux differentiable and its first derivative for £ £ ^"^(FjR 3 ) is 
given by 

dg[r,£]u = -[G(r)£,}G(r)u + (G(r)u ■ [0(r)£LV(r)) Af(r). 

Lemma 5.4. The mapping 

V: B°°(Q,s) -> ^(H s+1 (r),i7 s (r)) 
r i ^ r r divr r T r r 1 

is e 6" x ' -Gateaux differentiable and its first derivative for £ £ ^"^(TjR 3 ) is 
given by 

dD[r,£]u = -Tmce([G(r)$[G(r)u}) + ([g(r)u}M{r) ■ [G(r)£\Af{r)) . 

Similar results can now be obtained for the tangential vector curl by 
composition of the tangential gradient with the normal vector. 

Lemma 5.5. The mapping 

11: B°°{0,e) -> jSf(JJ s+1 (r),H s (r)) 
r i ^ T r curlp r t.,7 1 

is ff 00 -Gateaux differentiable and its first derivative for £ £ ^"^(TjR 3 ) is 
<jwen 6?/ 

<W2.[r, f]it = [g(r)t] T 7l(r)u - V{r)£ ■ TL{r)u. 

Proof. Let u e H S+1 (T). By definition, we have TZ(r)u = G{r)u A N(r). By 
lemmas 15.11 and 15.31 this application is ^^-Gateaux differentiable. For the 
derivative in the direction £ £ ^°°(r,R 3 ) we find 

d"fc[r,£]u = — [0 (r ) £] <? (r) uf\N(r) - Q{r)u A [0(r)ft/V(r). 

For any (3 x 3) matrix A and vectors b and c there holds 

(A6) A c + 6 A Ac = Tracc(A)(& Ac) — A T (b A c). 

We obtain the expression of the first derivative with the choice A = 
-[G(r)(i],b = G(r)uaiidc = Af(r). □ 

Lemma 5.6. The mapping 

TZ: B°°(0,e) -> J^(H s+1 (r), £F(r)) 
t* i ^ T r curlr r t,T 1 

is ff 00 -Gateaux differentiable and its first derivative for £ £ ^"^(FjR 3 ) is 
given by 

3 

d£[r, £]u = - ^ (£(r)& • tt(r)ui) - D(r)€ • ft(r)u 
i=i 

w/iere u = (wi,u 2 ,w 3 ) a«d £ = (6., 6, 6)- 
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Proof. Let ueH s+1 (r). Notice that we have curlr u = — Trace ([curlr u]j. 
We can therefore write 

n(r)u = -Trace(7l(r)u). 

The ^-differentiability of R results from the ^-differentiability of H. The 
first derivative in the direction £ is 

dH.[r, £]u = - Trace (dH[r, £]u) 

= - Trace ([£(r)£] T ["R.(r>]) - V(r)£ ■ Trace (-H(r)u]) 

3 

8=1 

□ 

Higher order derivatives of the tangential vector curl operator and of 
the surface scalar curl operator can be obtained by applying these results 
recursively. 

In view of the integral representations of the operators P^C^P" 1 and 
PrM^P" 1 , we have to study the Gateaux differentiability of the mappings 

r H> r r Ap 1 divr r t^ 1 
r H> r r Ap 1 curlr r TjT 1 . 

We have seen that for r g B°°(0,e) the operator curlr,. is linear and con- 
tinuous from H s+ (r r ) to H^(T r ), that the operator divr r is linear and con- 
tinuous from TH s+1 (r r ) to H^{T r ) and that Ap 1 is defined from H*(T r ) to 
H s+2 (T r )/C. To use the chain rules, it is necessary to prove that the deriva- 
tives at r = act between the spaces H s+1 (r) and H°(T) for the scalar curl 
operator, between the spaces TH s+1 (r) and H%(T) for the divergence opera- 
tor and between the spaces H*(T) and H S+2 (T)/C for the Laplace-Beltrami 
operator. An important observation is 

u r e H*(r r ) if and only if J r u r o (I + r) G H'(T). 

Using the duality (|2.5[) on the boundary T r we can write for any vector 
density j € H s+1 (r) and any scalar density <p £ H S (T) 

J T r (^cw\ Tr {Tr X j)} • <p J r ds = J curlr,. {t^j) ■ (r~ l (p)ds 

= J (r- x j) • curIr r (r-V) ds (5.1) 

j ■ r r f curlrv(T~V)) Jrds. 

Taking ip £ K (i.e. <p is a constant function) then the right-hand side vanishes. 
This means that J r (r r curlr r (t^j)) is of vanishing mean value. 
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Lemma 5.7. The mapping 

11*: B°°(0,e) -> ^(H s+1 (r),H:(r)) 
r i-> J r r r curlr r r^ 1 

is ff " -Gateaux differ entiable and we have in any direction £ = (£1,^2, £3) G 

^°°(r,R 3 ) 

or »=i 
<9 m 72* 

-g^rhe] = 0, for all m> 2. 

Proof. Looking at the expression of the derivatives of the tangential gradient 
and of the tangential vector curl in lemmas 15.31 and 15.51 we prove iteratively 
that all the derivatives of G(r){p and of H(r)<p are composed of Q(r)ip and 
TZ(r)ip, so that for ip £ K the derivatives of the right-hand side of (|5.1[) 
vanishes. We have for all m £ N and it £ H S+1 (F): 

— I ^ J r (r r curl rr (r-V)) ^} [r, £] = ^ ^ {ft* } [r, £]j = 0. 

It can also be obtained by directly deriving the expression of 72*ji using the 
formulas obtained in the lemmas 15 . 1 1 to [5~B1 The first derivative of 72.* is given 

by 

d (U*) [r, £]ti = -J(r) ■ E (0(r)& • 7t(r)u0 

i=l 

= -J r .r r ( E Vr^Cr" 1 ^) • curler" 1 !**) 

Vi=l 

The right-hand side is of vanishing mean value since the space Vr r H s (T r ) is 
orthogonal to curlrv H s (T r ) for the L (T r ) duality product. For the second 
order derivative we derive r H> d (72.*) [r, £]j in the direction 77 = (rji, 772, 773) <E 
^(r, R 3 ) and we obtain 



d 2 (ft*) [r; £, t}]ti = J(r) ■ E »)»»)& • TZ(r) Ul ) 

i=l 

3 



• E (Q(r)& • [G{rH T K{r 

i=l v 



= 0. 

Higher order derivatives of 72* vanish. □ 
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For the surface divergence, similar arguments can be applied. Using the 
duality (|2~3)l . we can write for j £ TH S+1 (T): 



T, 

r 



v(divrv(7r(r) 1 j)j ■ <p J r ds = j divr r (7r(r) 1 j) ■ (r r 1 ip) ds 



ds. 

This shows that for constant ip £ R the right-hand side vanishes and therefore 
J r (r r divr r (vr _1 (r)j')j is of vanishing mean value. We obtain the following 
result. 

Lemma 5.8. The mapping 

V* : B°°(0,s) -> ^(TH s+1 (r),HJ(r)) 

r 

is c £°° -Gateaux differentiable. 



i-» J r r r div Fr 7r 1 (r) 



Now it remains to analyze the inverse of the Laplace-Beltrami operator 
Ap. We apply the following abstract result on the Gateaux derivative of the 
inverse in a Banach algebra. We leave its proof to the reader. 

Lemma 5.9. Let U be an open subset of a Frechet space X and let y be a 
Banach algebra. Assume that f : U —¥ y is Gateaux differentiable at r$ £ U 
and that f(r) is invertible in y for all r £ U . Then g is Gateaux differentiable 
at ro and its first derivative in the direction £ £ X is 

df[r Q ^}=-f(r )- l odf[r ^}of(r Q )- 1 . (5.2) 

Moreover if f is c £ m -Gateaux differentiable then g is, too. 

From the preceding results we deduce the "if °°-Gateaux differentiability 
of the mapping 

C* : B°°(Q,s) -> jSf(iT s+2 (r),#!(r)) 

r i y JrTvArvr- 1 = -H*(r)H(r). 

Let us note that r r induces an isomorphism between the quotient spaces 
H s (T r )/C and H S (T)/C. 

Lemma 5.10. The mapping 

S°°(0, £ ) -> J?(H°(T),H s + 2 (r)/C) 

r h> (>C*(r)) _1 

is ff 00 -Gateaux differentiable and we have in any direction £ £ ff^iTjM. 3 ) 

d{r4(£'(r)) _1 }[0,(]= -Ap 1 o dC*[0, f] o Ap 1 . (5.3) 

Proof. We have seen in section [2] that the Laplace-Beltrami operator is in- 
vertible from H s+2 (T r )/C to H%(T r ). As a consequence C*(r) is invertible 
from H S+2 (T)/C to h{(T). We conclude by using Lemma [ESI □ 
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Let us give another formulation of (|5.3p . For any u £ H S+2 (T) and 
ip £ H~ S (T) we have 



r 



^Ar r (r,, 1 u)J ■ tpj r ds = — \Q{r)u ■ G{r)<pj J r ds 



It is more convenient to differentiate the right-hand side than the left hand 
side. For / £ H%(T r ), the element (£*(r)) 1 / is the solution u of 

-J (g(r)u-G(r)<p) J r ds = J f-ipds, for all ip £ H S (T). 

The formula (|5.3I) means that the first derivative at r = in the direction 
£ £ ^°°(r,R 3 ) of r i y (£*(r)) -1 / is the solution v of 

/ ^{( g ^ u °- g ^) Jr}[0,£]ds = - J Vr^-Vr^s, for all ip £ H S {T), 

(5.4) 

with uq — A r /. 

Now we have all the tools to establish the differentiability properties of 
the electromagnetic boundary integral operators and then of the solution to 
the dielectric scattering problem. 



6. Shape derivatives of the solution of the dielectric problem 

For the shape-dependent integral operators we now use the following simpli- 
fied notation 

(r) = * r E K P"\ * U „ (r) = 9 r Mlt P- 1 , 

C K (r) = P r C r K P-\ M K (r) = P r M r K P-\ 

In the following we use the results of the preceding paper [5] about the 
Gateaux differentiability of potentials and boundary integral operators with 
pseudo-homogeneous kernels. 

Theorem 6.1. The mappings 

B°°(0,e) -> if(TH-5(div r ,r),H(curl,X p )) 

r i-» V Mk (r) 

are infinitely Gateaux differentiable. The derivatives can be written in explicit 
form by differentiating the kernels of the operators ty r E and ty r M , see 
Theorem 4.7], and by using the formulas for the derivatives of the surface 
differential operators given in Section\^ The first derivatives at r — can be 
extended to bounded linear operators from TH2(divr,r) to H(curl,i7) and 
to H; oc (curl, VL C ). Given j £ TH 5 (divr,r), the potentials d^E K [0, £]j and 
d&M n [0,£,]j satisfy the Maxwell equations 

curl curl u — k 2 u = 

in Q and Q c , and the Silver- Miiller radiation condition. 
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Proof. Let j € TH 2 (divr, L) and let j — Vr£> + curlr q be its Helmholtz 
decomposition. Recall that ^ E K {r)j and ^M K { r )j can be written as: 

9 B „(r)j = ^ r -i( TrPr -ij) _ iv^r-^r.Ar^r-^)), 

* B „(r) j = curl^r- 1 ^,- 1 ^. 

By composition of differentiable mappings, we deduce that r i-> ^^(r) and 
r i— > $m s (0 are infinitely Gateaux differentiable far from the boundary 
and that their first derivatives are continuous from TH 2 (divr, T) to L 2 (17) U 
L loc (ri c ). Recall that we have 

curl* E K (r)j = n^M K (r)j and curl ^M K (r)j = n^ Eii {r)j. 

Far from the boundary we can invert the differentiation with respect to x 
and the derivation with respect to r, which gives 

curl d^E K [0, £]j = k d$>u K [0, £] j and curl d* M(i [0, £] j = « d* [0, 

It follows that d*E„[0,£]j and d*Af„[0,^]j are in H(curl, Q) U H loc (curl, IF) 
and that they satisfy the Maxwell equations and the Silver-Muller condition. 

□ 

We recall from Section |4] that with the notation of Section[5]the operator 
C K (r) admits the following representation 

C K (r)j =P r C r K P r - 1 j = V r F(r) + curl r Q(r), (6.1) 

where 

P(r) = -^rWj-^'W^KV 1 ) [5(r)p + 7l(r)g] 

and 

Q(r) = k {C*{r))- x V*(r)<K(r)T- x (tvVJt-- 1 ) [Q(r)p + H(r)q] 
+l{ Tr V^T- 1 ) (r r A Vr {T- l p)) . 

r\j 

Let j £ TH s (divr,r) and let j = Vr p + curlr q be its Helmholtz 
decomposition. We want to derive 

PrC^j = P r C£(Vr r T-V + curlr,. r^q) 
= P r (Vr r P r + curl Fr Q r ) 
= Vr-P(r) + curlr Q(r). 

We find 

dC K [0^}j = V r dP[0,e] + curl r dQ[0,$. 

Thus the derivative with respect to r of P r C^P r _1 j is given by the derivatives 
of the functions P{r) — r r {P r ) and of Q{r) = r r {Q r ). 

We also note that for an r-dependent vector function f(r) on T there holds 

rf{ 7 r(r)r- 1 /(r)}[0,e]=7r(0)d/[0^]. 

By composition of infinitely differentiable mappings we obtain the fol- 
lowing theorem. 
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Theorem 6.2. The mapping: 

B°°{0,e) -> jgf (TH^^(div r ,r),TH-5(div r ,r)) 

r i — ^ PrC^P" 1 

is infinitely Gateaux differentiable. The derivatives can be written in explicit 
form by differentiating the kernel of the operator C r K , see Corollary 4.5], 
and by using the formulas for the derivatives of the surface differential oper- 
ators given in Section^ 

Similarly, recall that the operator P r M£P r -1 admits the following rep- 
resentation : 

PrM^Pr 1 ] = Vr-P'(r) + curlr Q'(r), 

where 

P'(r) = (C*(r))- l (K*J r -T r nr ■ (tv^t" 1 ) [Q(r)p + H(r)q] 
+(£*(r))- 1 (J r • T r D:T-i)(T r A rr (T-ip)) 

and 

Q'(r) = {C*{r))^n*{r){r r {Bl - Z^r" 1 ) [G(r)p + n{r)q] 

with 

TrBlP^j = r r y VG{k, I • -y r \) (n r ( • ) • (Vr^- 1 ^)) 

+ jf VG(k, I • -ifrl) (n r ( ■ ) ■ (curlr r Ty^Xifr)) ds(»r)}| . 

Theorem 6.3. The mapping: 

B°°{0,e) -> ^f(TH-^(divr,r),TH-5(div r ,r)) 

«s infinitely Gateaux differentiable. The Gateaux derivatives have the same 
regularity as M K , so that they are compact operators in TH~ 5 (divr, T). The 
derivatives can be written in explicit form by differentiating the kernel of 
the operators M£, see (HI Corollary 4.5], and by using the formulas for the 
derivatives of the surface differential operators given in Section [21 

Proof. The differentiability of the double layer boundary integral operator 
is established in Example 4.10]. It remains to prove the infinite Gateaux 
differentiability of the mapping 

B s -> if (TH-5(div r ,r),H5(r)) 
r ^ r r B^p-\ 

The function (x,y — x) M> VG(k, \x — y\) is pseudo-homogeneous of class 0. 
We then have to prove that for any fixed (x,y) € (r x T)' and any function 
p G (r) the Gateaux derivatives of 

r i— > (T r n r )(x) ■ (r r 'Vr 7 .T~ 1 p) {v) 
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behave as \x — y\ 2 when x — y tends to zero. To do so, either we write 

(r r n r )(x) • (-rvVrvr^V) (v) = ((r r n r )(x) - (r r n r )(j/)) • (r r Vr r r r "V) (y) 
or we use Lemmas 15.11 and 15.31 and straighforward computations. □ 

Theorem 6.4. Assume that E mc £ H 1 1 oc (curl, R 3 ) and that the mappings 

B°°(0,e) -> ThH(div r ,r) 
r P r (n r A E|r°) 

r h> P r (n r A (curlE mc ) |r J 

are Gateaux differentiable at r — 0. TTien i/ie mapping from r £ B°°(0, e) to 
the solution $(r) — E(f2 r ) £ H(curl, f2 r ) U Hi oc (curl, f2 c ) o/ ifte scattering 
problem by the obstacle £l r is Gateaux differentiable at r = 0. 

Proof. We use the integral equation method described in Theorem 13.11 Let 
j be the solution of the integral equation p.8j) . By composition of infinitely 
differentiable mappings we see that 

B°°(0,e) -> J^(TH-^(div r ,r),TH-5(div r ,r)) 
r h> S(r) = P r S r p- 1 

is "if °°-Gateaux differentiable. Then with (|3.6[) and (|3.8p we get for the exte- 
rior field S s 

dS s [0, £] = (-dttj^ [0, £] - «^*Af Ke [0, £]C * - dC * [0, ei) j 

+ (-*e Kc - iv^M^C^S- 1 ^- dS[0,£] j) 
+ i^M^Co^S- 1 (-^M K j0,C]7i?E mc - daj(U]7^E mc ) 

+ (-3^ - i^.CSJS" 1 + M Ki ) d{P r7 LE mc } [<U]) 

+ (-* Bks - i^^CJJS- 1 (-C^rfjp^E^} [0,£]) . 

We know that j £ TH2(divr,r), so that the first terms on the right-hand 
side are in Hi oc (curl, f2 c ), and the hypotheses guarantee that the last two 
terms are in Hi oc (curl, fl c ). For the interior field we write 

d£% £] = - i«f* B [0, e]7k (E* + E"") - rf*M». [0, Hhh (E s + E mc ) 

--*^d{ P ^. (^M + E mc )}[0,^] 

-Vm^A^d {<? s (r) + E inc )}[0,Z}. 

The hypotheses guarantee that 7^ (E s + E™ c ) and 7^ (E s + E inc ) are in 

TH2(divr,r), which implies that the first two terms are in H(curl, f2), and 
that the last two terms are in H (curl, Q) . □ 
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Theorem 6.5. The mapping sending r G B°°(0,e) to the far field pattern 
E°°(51 r ) £ T&°°(Sr) of the solution to the scattering problem by the obstacle 
fl r is -Gateaux differentiable. 

Proof The mapping B°°(0,e) 3 r h-> {(x,y) h-> e lRi <y +r ^} G c if°°{S 2 x 
r) is ^"^-Gateaux differentiable and the derivatives define smooth kernels. 
By the linearity of the integral we deduce that the boundary-to-far-field 
operators 

B s -> J?{TH s (dW r ,T),tf°°(S 2 )) 



are ^ °°-Gateaux differentiable. For j G TH s (divr, T) we have: 

dVgJO,Z]j(x) = i K xA ^e- iK *y(div r t;(y)-iKX-Z(y))j(y)ds(y?) Ax, 
and 

rf^/jO,C]j(x) =kxA (^e- iK *y(div r Z(y)-iKx-Z(y))j(y)ds(y)) . 

We conclude by using the integral representation of E°°(f2 r ) and previous 
theorems. □ 

6.1. Characterization of the first derivative 

The following theorem gives a caracterization of the first Gateaux derivative 
of r i— !► S(r) at r = 0. 

Theorem 6.6. Under the hypotheses of Theorem 16.41 the first derivative of 
the solution <o{r) of the dielectric scattering problem at r = in the direction 
£ G ^""(FjR 3 ) solves the following transmission problem : 

( curl curl dS l [0, £1 - n 2 d§ 1 [0, £1 = 

<^ I (6.2) 

[ curl curl d<? s [0,£] - K 2 e d£ a [Q,£\ = 

with the interface conditions 

( n AAT[0,£] -nAd<? s [0,£] = g D 

\ A curl d<T [(),£] - (i- 1 nAcwrldg s [0,£\ = g N , 
where with the solution (E 1 , E s ) of the scattering problem, 

g D = - (£ • n) (n A curl E'-nA curl(E s + E rnc )j A n 

+ curlr ((£ • n)(n • E* - n ■ {E s + E mc )f 

and 



g N = — (£ • n) ( —n AE' -n A (E s + E mc ) ) A n 



+ curlr ((£ • n) curl r E 4 - M J X curl r (E s + E mc )) 

and rf(? s [0,£] satisfies the Silver-Muller radiation condition. 
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Proof. We have shown in the previous paragraph that the potential operators 
and their Gateaux derivatives satisfy the Maxwell equations and the Silver- 
Miillcr radiation condition. It remains to compute the boundary conditions. 
They can be obtained from the integral representation, but this is rather 
tedious. A simpler way consists in deriving for a fixed x € T the expression 

n r {x + r{x)) A (<T(r)0 + K*)) ~ + r ( x )) - E " lC ( x + r ( x ))) = 0. 

(6.4) 

This gives 

= (W[0,£](x) A (E*(a;) - E s {x) - E mc {x)) 

+ n(x) A {dS l [0, i\ (x) - dS s [0, £] (x)) 

+ n A (£(x) • V (E* - E s - E mc )) . 

We now use the explicit form of the shape derivatives of the normal vector 
given in Lemma [5.1 1 : dftf[0, £] = — [Vr£] n, and the formula Vu = Vr« + 
(f£) n. We obtain 

n(x) A (d<T* [0, £] (a:) - di a [0, f ] (x)) = 

[V r e]nA(E l ( a; )-E s (x)-E mc ( a ;)) 
- n A • V r (E*(x) - E s (z) - E^fi))) 

- (C • n)n A A (E\x) - E s (x) - E mc {x)) . 

Since the tangential component of E l — E s — E mc vanishes, we have 

(Z{x)-V T {E\x)-E s {x)-E mc {x))) 

= (n-(E l (x) - E s (x) - E mc {x))) ([v r n T ] £ 

and 

(lV r $n)A(E\x)-E s (x)-E mc (x)) 

= ([V r f] n) A n (E 4 (x) - E s (x) - E inc (a;)) • n 

Since we are on a regular surface, we have Vr« = Vr^ T and 

([Vrf] i) A n - n A ( Vr« T £) = curl r (£ ■ n) . 



Using the expansion (see [HI p. 75]) 

curl m = (curlr «) n + curlr («r • n) — ([Vrwlu) An— — — 

\an 

we obtain that 

n A ( 7n E* - 7^(E S + E mc )) = — n A (7 curl E* - 7 c curl (E s + E mc )) A n 

+ curlr (n- ( 7 E i - 7 c (E s + E inc ))). 
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Here we used that the curvature operator [Vr**] acts on the tangential com- 
ponent of vector fields, so that 

[V r n] ( 7 E i - 7 c (E s + E inc )) =0. 

Thus we have 

9 D =-(f.n)nA (7 curl E 4 - 7 C curl (E s + E mc )) A n 
+ curlr ((£ • n) (n ■ 7 E l - n • 7 C (E S + E" lc ))) . 
To obtain the second transmission condition, we use similar computations 
with the electric field E replaced by the magnetic field curl E. □ 

7. Perspectives: Non-smooth boundaries 

We have presented a complete differentiability analysis of the electromagnetic 
integral operators with respect to smooth deformations of a smooth bound- 
ary in the framework of Sobolev spaces. Using the boundary integral equation 
approach we have established that the far-field pattern of the dielectric scat- 
tering problem is infinitely differentiable with respect to the deformations 
and we gave a characterization of the first derivative as the far-field pattern 
of a new transmission problem. 

In the case of a non-smooth boundary — a polyhedral or more generally 
a Lipschitz boundary — the formulas determining the first derivative given 
in Theorem l6.6l are problematic. The normal vector field n will have discon- 
tinuities, and the factor £ • n and vector product with n that appear in the 
right-hand side of (|6.3|) may not be well defined in the energy trace spaces. It 
is, however, known for the acoustic case that the far field is infinitely shape 
differentiable for non-smooth boundaries, too, see |19) for a proof via the 
implicit function theorem. 

Our procedure of using a boundary integral representation gives an al- 
ternative way of characterizing the shape derivatives of the solution of the 
dielectric scattering problem and of its far field. We do not require the compu- 
tation of the boundary traces of the solution and taking tangential derivatives 
and multiplication by possibly discontinuous factors. Instead we determine 
the shape derivatives of the boundary integral operators. While the study 
of Gateaux differentiability of boundary integral operators, for the case of 
smooth deformations of a Lipschitz domain, is still an open problem that 
will require further work, our approach via Helmholtz decompositions seems 
to be a promising starting point for tackling this question. Let us briefly 
indicate why we think this is so. 

We consider the case where the boundary F is merely Lipschitz, but 
the deformation is defined by a vector field £ that is smooth (at least < ^' 1 ) 
in a neighborhood of F. Note that the reduction to purely normal displace- 
ments that is often used for studying shape optimization problems for smooth 
boundaries does not make sense here, as soon as there are corners present. 
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In this situation, many of the ingredients of our toolbox are still available. 
Here are some of them. 

First, the change of variables mapping r r still defines an isomorphism be- 
tween H 2 (r r ) and H* (r). By duality, we see that the mapping u r i— > J r r r u r 
defines an isomorphism between H~z{T r ) and H~^(T). When we want to 
transport the energy space TH~ 5 (divr r , r r ), we can still use Helmholtz de- 
composition. Namely, the following result is known, see [U EJ El HI [5] . 

Lemma 7.1. Assume that Y is a simply connected closed Lipschitz surface. 
_ i 

The Hilbert space THn 2 (divr,T) admits the following Helmholtz decomposi- 
tion: 

THj[*(div r ,r) =V r K(r)fficurl r ili(r). 

where 

H(T) = {ueH 1 {T) : Aru G H~^ (Y)}. 

_ i 

The notation THn 2 (divr, r) recalls the fact that special care has to be 
taken for the definition of the energy space. 

A natural idea for the transport of the energy trace space is then, instead 
of (ETo) . to define 

P r : TH~*(divr r ,r r ) — > TH^(div r ,r) 

Vr r Pr + curlr,. q r i-> Vr Ap 1 (j r (r r Ar r p r )) + curlr(T r g r ). 
This is justified by the sequence of isomorphisms 

H(T r )/C ^ H**(T r ) — > H^(T) U(T)/C 

Pr !-> Ar r p r !-> Jr(TrAr r Pr) ^ Ap 1 ( J r (T r Ar r Pr)) ■ 

The inverse of the transformation P r is given by 

P.- 1 : THp(div r ,r) — ► TH~*(divr r ,r r ) 

Vrp + curl r q i ^ V r ,. r r T 1 (£*(r)) 1 A r p + curl r ,. (r^ 1 ^). 

In this situation it seems to be more convenient to rewrite the operators 

PrC^PjT 1 and PrM^P^ 1 as operators acting on the the scalar fields p* — 

_i i 
Ar\p G -ff* 2 (r) and g G H^(T r )/C instead of p and q. For example, the 

operator P r C r R P r is defined for j = VrA r p* + curlr q G TH | 2 (divr,T) 

by 

P^P; 1 = VrA? J P*(r) + curl r Q(r), 

with 

P(r) = -/c^W^r- 1 ) [e(r)(£*(r))"V+^(r)g 

and 

Q(r) = -« (£*(r))- x X>*(r)7r(r) (r^T" 1 ) [0(r)(£*(r))~V + tt(r)g 
+ \ (TrVZr-i) (J-V) . 

fx 
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Here we have used the same notation for the surface differential operators as 
introduced in Section [5] These formulas together with similar ones for the 
operator P r .M re P~ can now be the starting point for generalization of the 
analysis of shape differentiability of the Maxwell boundary integral operators 
to Lipschitz domains. We expect that the results for the differentiability of the 
surface differential operators and then also of the boundary integral operators 
will be similar to what we have obtained for the case of smooth domains. This 
is, however, far from trivial and will require further work. 
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